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RESEARCH MEMORANDUM 


BOUNDARIES OF SUPERSONIC AXISYMMETRIC FREE JETS 

By Eugene S. Love, Ifl.ldred J. Wood ling, 
and Louise P. Lee 


SUMMARY 


Calculations have "been made "by the method of characteristics of 2960 
boundaries of supersonic axlsymmetrlc free jets esdiaustlng frcm conically 
divergent nozzles Into still air. Ihe calculations cover jet Mach num- 
bers from 1.5 to 3 * 0 j ratio of specific heats of the jet from I.II5 
to 1.667, divergence angles of the nozzle from 5*^ to 20°, and jet static- 
pressure ratios from 1 to 10 . Hie results are presented nondlmenslonally 
and In a form that facilitates practical lOse. Some random examples of 
the effects of the variables Involved are shewn, and the possibilities 
of the simulation of one jet by another are examined for jets exhausting 
Into still air and Into si:^ersonlc streams. 


INTRODUCTION 


One of the important ctiaracteristlcs of free jets exhausting into 
still air is the shape of the jet boundary. Once the shape of the bound- 
ary is known. It may be applied to a number of problems for which the 
condition of constant pressure along the jet boundary Is a permissible 
as8un5)tlon. Some of these applications have been dlsciossed and Illus- 
trated in reference 1 . Perhaps the most direct applications occur in 
the design of airframe conponents to obtain jet clearance in the static 
thnist condition and In the assessment of the deflections of the external 
stream caused by a propulsive jet exhaiistlng from an airplane flying at 
low to moderate subsonic speeds. Other applications occur In estimating 
the performance of ejectors having near-constant pressure xlng and l-n 
the estimation of pressure on the base annulus separating an exhausting 
jet from an external sipersonlc stream. (See ref. 1 .) 

In reference 1 several calculations were made by the method of 
characteristics to obtain primarily the shape of the jet boundary, and 
It was shown that the theoretical nonviscous botindarles agreed satis- 
factorily irLth observed experimental boundaries over the region of prin- 
cipal interest, namely, from the nozzle exit to the vicinity of the 
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maxl.TmTm diameter of the free Jet within the first wavelength. Beyond the 
maxiimnn diameter, turbulent diffusion significantly affects the boundary 
shape and generally permits simplified assun^itions to he made of the shape. 
Hhe theoretljcal boundaries of reference 1 included a few calculations for 
‘conically divergent nozzles in which the effects of nozzle divergence 
angle and ratio of specific heats of the jet were examined for certain 
arbitrarily chosen initial conditions. Ghe purpose of this paper is 
to extend the theoretical calculations of boundary shapes for jets 
exhaiisting from conically divergent nozzles and to present these bound- 
ary shapes in a form convenient for practical use. Ghese boundaries 
cover the following range of variables: jet Mach numbers from I.5 to ^. 0 , 

jet static-pressure ratios from 1 to 10, nozzle divergence angles from 
5° to 20°, and ratio of specific heats of the jet from I.II5 to I.667. 


SYMBOLS 


See figure 1 for clarification of symbols. 


M, 


jet Mach number; value at nozzle lip immediately ahead of exit 


Poo 


jet static pressure; value at nozzle lip immediately ahead of 
exit 

ambient static pressure 
radius of jet exit 


X distance from plane of nozzle exit measured parall el to 

jet axis 

y peipendicular distance frcan jet axis 


ratio of specific heats of jet 

initial inclination of free jet boundary at jet exit 


% 


conicELl divergence angle or exit angle of nozzle measxired with 
respect to jet axis 


Subscripts : 

5 jet 


N 


nozzle 
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o value at 0jf = 0° 

00 anibieiit 


CALCULATIONS 


Scope 

Q3ae calculations were made for all combinations of tLe following 
varlaLles: values of Mj of 1.5^ 2.0, 2.^, and 5*0j values of /j 

of 1 . 115 , 1.200, 1 . 300 , 1.400, and 1.66'Jj values of Bjj of 5 °, 10 °, 
15°j and 20°j and values of of 1 to 10 in increments of O. 25 . 


Msthod 

The method of characteristics was used throughout. Bie procedure 
TOs identical to that described in the section of appendix A of ref- 
erence 1 entitled "Solution for Divergent Nozzles." 

The leading characteristic lines (fig. l) and the initial values 
for the characteristic lines from the center of expansion were computed 
manually. Bae characteristic nets were determined on a card-programed 
electronic calciilator. 

ihe calculations were proposed to he carried throu^ only until 
the maxi tmitn diameter of the Jet boundary was reached, since the shape 
beyond the maxj.imim diameter has little practical value. It will be 
noted, hOTirever, that a nmber of the boundary calculations extend beyond 
the maxiinum diameter. Biis random extension results from a random con- 
tinuance of the calciilations; since these additional points on the 
boundary were available, they have been included. For those boundaries 
that extend well beyond the maximum diameter, the boundary shape in the 
nei^borhood of the final calciilated point on the boundary may be sub- 
ject to small errors that the foldback method of cdiaracteristlcs c^ 
introduce in this region. (See ajpendix A of ref. 1.) A few of the 
boundary calculations at Mj = I .5 were inadvertently tenninated slightly 
short of the maximum diameter. . In all the present calculations the first- 
family characteristic line determining the final point on the boundary 
had its origin on the leading-characteristic line (fig. l) , Consequently, 
the leading-characteristic line was terminated short of the Jet axis, and 
the conplications of calculations ccmlng off the axis were not encount- 
ered. (See appendix A of ref. 1.) Inasmuch as the accuracy of the 
boundary calculations is closely related to the density of the char- 
acteristic net, some idea of the number of points contained in a typical 
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characteristic net (see fig. 1 ) may he of interest. For every increase 
in — of 0.25, an arLditional characteristic line from, the center of 


expansion ■was in-fcroduced. Bans, a relati-vely uniform accuracy "was main- 

Pi o 

tained as ^ "was increased. For exanple, at = 1 * 5 ^ Qjj = 15 .> 


and = l.li-OO, there are 28 points in the calculation at 

p.. 

and 288 points in the calcu l ation at — = 10. 




= 1 


Qhe calculated houndaries are presented in discontinuous-slope form, 
which is conpatihle -wl-th the la'tbice-point characteristic method (that 
is, calcvilated points on houndary Joined hy s-trai^t lines) . In practical 
application a smoo-th cirrve is, of coirrse, fit-ted to the discontinuoi:is- 
slope form of the boundary, particular care being given to maintain -the 
initial inclination of the ho'undary. 


RESULTS 


Presentation of Boundaries 


Bie 2960 calculated hoxiudaries are presented in figures 2 to 5. 

Each group of houndaries corresponds to a -variation of ^ from 1 to 10 


for a particular combination of Mj, 9 jf, and /j. Ihe solid ho\mdaries 


deno-te changes in 


El 


of 1; -the three dashed houndaries between the 


pj 


solid houndaries correspond to changes in 


of 0.25 and are elIso cal- 


culated not interpolated houndaries. For clarity, a number of the groups 

ha've insets covering -tiie lower range of — in which the ordinate ^ 

is shown to -twice scale. The order of presen-fcation of -these groups of 
ho^lndaries is as follows; 


Figure 


2 

1.5 

3 

2.0 


2-5 

5 

5.0 
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•For figures 2 to parts (a), (t>), (c), and (d) are for 0jj 
of 5°^ 10°, 15°^ 20°, respectively. Each figure subdivision contains 

five groups of houndaries corresponding to the five values of /j, and 

these five groups are presented in the order of ascending values of 7^, 
that is, fran 7j = l.U-5 to I.667. 


Effects of Varlahles 


Qhe effect of increasing static-pressure ratio — is ohvlous 

^00 

from figures 2 to 5* Effects of the other variables at a static-pressure 
ratio of 5 (arbitrarily chosen) are shovm in figure 6. At this static- 
pressure ratio, the effect of Increasing Mj (fig:. 6(a)) is to decrease 
the initial Inclination of the boundary, to increase the maxi.muni diameter 
of the free jet, and to move the max^mimi diameter farther away from the 
plane of the jet exit. The effect of increasing 7^ (fig. 6(b)) is to 
decrease the initial inclination of the boundary, to decrease the max- 
imum diameter of the free jet, and to move the maximum diameter closer 
to the plane of the jet exit. Ihe effect of increasing 9jj (fig. 6 (c)) 


is to increase the initial inclination of the boundary, to increase the 
mR-x-tTm-RTi diameter of the free jet, and to move the maximum diameter closer 
to the plane of the jet exit. A brief examination at other values 


of 



indicated that these effects are apparently typical of all 


pressure ratios; however, there ml^t be exceptions which were not dis- 
closed by this examination. 


Simulation of One Jet Boundary by Another 


In figures 2 to 6 it may be noted that the initial inclination of 
the boundary varies with every variable. When jet boundaries are simu- 
lated, the initial inclination is, perhaps, the most Important property 
to be duplicated. (See refs. 1 and 2.) For this reason and becatise of 
the problems confronted in jet simulation, it is of Interest to examine 
the degree of boundary simulation that may be achieved by simulating 

P-? 

initial inclination only throu^ a variation of — only and, in addl- 

^eo 

tion, to determine whether conplete simulation is possible throu^ a 
variation of all variables. 


Pj 

In order to aid in the selection of values of — that give the 

CO 

same initial inclination 6j, figure 7 has been prepared which gives the 
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P-I 

variation with — of the initial inclination for a Jet exhausting from 

^00 

a nozzle having zero divergence angle ^denoted hy . Qlie value 

of at any finite value of % is s i m ply ' the sum of 0]j + 

Figure 8 presents some typical examples of boundary simulation. 
Kandom examples of the hoiindary simulation that may he ej^ected hy 
duplication of only are shown in figures 8(a) to 8(d) . In these 

J Pi 

examples and 9 jj remain consteuat and is allowed to vary so 

■^00 

that the value of 8j is "^e same for hath values of 7^. (See fig. J.) 

In essence, a houndary given hy an air jet (7^ = 1 . 400 , dashed houndar;^ 

is compared with a houndary given hy a jet of another gas or hy a heated 
air jet (7^ < 1 ^ 0 (^ . Hie resiilts indicate that the differences between 

the boundaries thus simulated can he attributed primarily to differences 
in 7j and that differences in 7^ of the order of 0.1 or less create 
negligibly small differences in the boundary when the Initial inclination 
is duplicated. Th i s observation also implies that, insofar as the effects 
of hoTindary shape are concerned, air jets heated to the order of 2,500° R 
or lower can he satisfactorily simulated hy cold air jets hy altering 
the pressure ratio to give the same initial inclination. Hie results 
discussed in the appendix show that this observation may he extended to 
the case of supersonic jets exhausting into supersonic streams. 


Figure 8(e) gives an example which illustrates the point that the 
boundary of any jet may he satisfactorily simulated hy an air jet hy 
allowing a freedom of choice of the variables for the air jet. Extremely 
Pi 

large values of — would, of course, produce exceptions. (It will he 


noted that the boundaries of figure 8 are, as before, given in the 
discontinuous-slope form. Hiis form tends to exaggerate the differences 


in the lower range of 


— , particularly at the first calculated point on 


the boundary.) 


Figure 8(f) illustrates that, for a particular value of 7^, there 

Pi 

are a number of comb j nations of Mj, 0 jj, and — that produce essen- 
tially the same boundary. 


Interpolation and Extrapolation of Calculated Boundaries 


Interpolation and cross -plotting can he used to obtain a reliable 
boundary ifithin the range of any of the variables covered by these cal- 
culations . 
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The reliability of a boundary necessitating the use of extrapolation to 
cover a varialile the magnitude of ■which places it outside the range of 
■fchese calculations is not readily apparent. A cursory examination, ho\r- 


Pi 

ever, gave the following indications. Extrapolations of — should not 

p. ^00 

exceed -values of — of about llf and probably less; extrapolations 


p 

below — = 1 are 'not recommended. Extrapolations to 0jj = 0° are of 
^00 


particular Interest since the nozzle having parallel flow at the exit 
is often encountered. This extrapolation may be made -with good results 
as Indicated by the' exan^jle given in figure 9 vhere an extrapolated 
boundary for 0j^ = 0° is shown to be in good agreement with the cal- 
culated boundary for 0jj = 0° from reference 1. Extrapolations to 
v^ues of in excess of about 25° are not recommended. Extrap- 
olations of Mj to Mj = 1 (% = 0°) appear permissible, and extrap- 


olations to Mj 5.5 or sli^-biy hi^er appear within reason, ihe 
necessi-ty for extrapolations of • beyond the range of these calcula- 
tions appears remote; however, should the need arise, ilie extrapolation 
could in all likelihood be made 'wl-th confidence. 


CONCLUDING REMARKS 


Calciilatlons ha^ve been mside "by -the method of characteristics of 2960 
boundaries of si: 5 )ersonlc axlsymmetric free jets exhausting from conically 
divergent nozzles into still air. The calculations cover jet Mach num- 
bers from 1.5 to 5*0^ ratio of specific heats of the jet from I.II5 
to 1 . 667 , divergence angles of the nozzle from 5 ° 'to 20°, and jet static- 
pressinre ratios from 1 to 10. Ihe results are presented nondlmensionally 
and in a form -that facilitates practical use. Scane random exangJles of 
the effects of the variables involved are shown, and the possibilities 
of the simulation of one jet by another are examined for jets exhausting 
into still air and into sig)ersonic streams. 


Langley Aeronautical Laboratory, 

Natlonai Ad-vlsory Committee for Aeronautics, 
Langley Field, 1b.,, June 22, 1956. 
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APPENDIX 

JETS E5fflAUSTING INTO SUPERSCaHC STREAMS 


Because of the negllglhly small differences found, in the houndaries 
of supersonic Jets esshausting into still air for differences in 7j of 

the ord.er of 0.1, provided the initial inclination is duplicated, an 
attenpt was made to determine whether equally small differences in the 
houndary shape and, of equal Inportance, in the ambient wave inter- 
ference field could he expected for stipersonic Jets exhaiistlng into super- 
sonic streams under the same conditions of simulation (identical Sj). 


At the outset, one mi ght suspect that the effect of a given change 
in would have less effect upon the Jet-houndary shape for Jets 

exhausting expansively into supersonic streams as conpared with Jets 


ejdiausting into still air at the same value of 
for the following reasosas. For the same value 




Blis suspicion arises 
, the Jet exhausting 


into the supersonic stream undergoes considerahly less eijqpanslon than the 
Jet exhausting into still air because of the pressure rise across the 
exit shock generated in the supersonic stream. Thus, effectively, the 
initial portion of the issioing Jet flow for the Jet exhausting into the 
supersonic stream resembles that which occurs for a Jet exhausting into 

p ' 

p. 

has already been shown that the lo\fer the value of the less the 

Pi 

7j, and that at ^ = 1 changes in 7j greater than 0.1^• 


still air at a significantly lower value of 


In reference 1, it 


effect of 


have negligible effect upon the boundary shape of Jets exhausting into 
still air. 


A characteristic celculatlon has been made by Schafer in reference 5 
for a Jet having = 2.57^ 7j = 1-^05, = 12. 5°^ and exhausting 

P1 

at — » 8.9 into an ambient stream having a Mach number of 5.2k and 


a ratio of specific heats of l.k05 (same as Je'fc); "the exterior surface 
ahead of the Jet exit had a boattall angle of reference 4, the 

isobars for this calculation have been ccsiputed and the calculation is 
presented in isobar-streamline form. (The values of some of the initial 
variables quoted in figure I5 of reference 4 contain slight errors which 
were corrected for this calculation.) The Schafer calculation has been 
repeated with the initial variables remaining the same except for 7^, 
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Pi 

which was changed to 1.300^ and which was reduced so that the InltlaJ. 

^00 

lncl±nation of the mixing boundary was duplicated. This new calculation 
was carried to the TnaxiTmim diameter of the jet. The houndary shape and 
the interference Isohars in the ambient stream obtained in this cal- 
ciilation could not be distinguished from those shown in reference 4. 

Thiis, the conclusion dra™ for jets exhaiosting into still air would appear 
to apply as well to jets exhausting into svipersonic streams, namely, 
differences in 7^ of the order of at least 0.1 create negligibly small 

differences when the initial inclination of the jet boundary is duplicated. 

On the basis of the above conclusion, cold air jets may be expected 
to simulate the wave interference of other jets satisfactorily in the 
range of variables covered in reference 2 when the initial inclinp.tion 
is duplicated according to the curves of reference 2. (The restilts of 
reference 2 are for /j = 9/T or 1.286; therefore, the cold-air- jet 

simulation amounts to a difference in 7^ of the order of 0.1.) 

Furthermore, the curves of reference 1 giving Initial inclination of the 
boundary and of the exit shock may be used to achieve simulation of the 
wave-interference field with cold air jets over a larger region of the 
interference field than was at fihst thou^t possible. 
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Final calculated point on boundary 
Region covered by characteristic net —j 



Figure 1,- Sketch of nozzle. Jet hoimdary, and pertinent Information^. 




Flgiire 2.- Jet boundarleB at Mj = I .5 for Jet static-pressure ratios 
from 1 to 10. (Dashed boundaries denote changes in of 0.25.) 
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Figure 2.- Contlmied. 
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(a) 9 ]j a 5 °. ContlJiiied. 


Plgur© 2.- Continued. 
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(t) 0JJ = 10°. Continued.- 


Figure 2.- Continued 
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("b) 0N “ 10°. Continued. 
Figure 2.- Continued. 








(b) 0 ]j a l(r . Continued . 

Figure 2 .- Continued. 
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(b) 9jj a lo' 


Figure 2 



X 

n 


Concluded. 


Continued 


[3 












(c) 0 uf “ 15°* Continued. 


Figure 2 .- Continued 
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Figure 2.- Continued. 






Figure 2.- Continued. 




(c) Sfi = 15° • Ctoncluded. 


Figure 2 .- Continued 
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(d) 6j^ »» 20°. Contimied. 

Figure 2.- Continued. 
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(d) 0JJ = 2CP . Concluded. 
Figure 2.- Concluded. 
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(a) 0N = 5°. 

Figure 5-” Jet bo undar ies at Mj =2.0 for Jet static-pressure ratios 
from 1 to 10. (Dasbed boundaries denote changes In Pj/Poo 0*25.) 
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(a) Bjj a 5°. Continued. 
Figure 5*- Continued, 
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(a) 0 u a 5 °* Contltitisd. 


Figure J.- Continued 
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(a) 0ijf “ 5°- Concluded. 


Figure 3.- Continued. 
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Figure 5.- Contlnujed . 
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(Td) 9 if = 10 °. Continued. 
Figure 3.- Continued. 
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("b) 0JJ “ 10°. Continued. 
Figure 5.- Continued. 
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(b) 0|j 10°. Continued. 

Figure 5*~ Continued. 
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Figure 5.- Continued. 
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(c) 0-xT = 15° • Contlniied. 
Figure 5 .- Continued, 
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Figure 5*“ Continued. 
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(c) 6 jj = 15°. Continued. 
Figure 5 -“ Continued. 










(d) 0JJ 20^ • 

Figure 5.- Continued. 
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(d) 0JJ = 20°. Continued , 
Figure 5 .- Continued. 
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Figure 5"~ Continued 
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Figure 4.- Jet ‘bcrurtdarles at Mj = 2.5 for Jet static -presstire ratios 
from 1 to 10. (Dashed 'boundaries denote changes in Pj/Po, of O.25.) 
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(a) 9 jj o Continued. 

Figure 4.- Continued. 







(a) Qjj = 5°» Continued. 
Figure il-.- Continued. 
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(a) Bjj = 5°. Conclvided. 
Figure li-.- Continued. 
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(ij) 0JJ ■= 10*^. 

Figure 4 .- Continued. 
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(d) Qu ~ 10°' CorLtinued . 
Figure 1;,- Continued. 
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(c) 0 JJ = 15 °. 
Figure Ij-.- Contlnujed. 
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(c) 0N = 15 ° • Continued. 
THnnire Continued. 
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(d) 0JJ = 2CP. 
Figure 4.“ Continued 




Figure Ij-.- Continued 
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Figure if.- Contla\ied 
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(d) 0J]- = 20°. Continued. 

Figure k.- Continued. 
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Figure 5*- Jet bou n daries at Mj - 5*0 for Jet static -pres sure ratios 
from 1 to 10. (Daalied boundaries denote changes in Pj/Poo 0 * 25 .) 
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Figure 5*“ Contimed. 
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("b) 0jf = 10°. Continued 
Figure 5 .- Continued. 
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Figure 5 *~ Continued. 
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(c) 0jf = 15°* Continued. 
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Figure 5*’’ Continued 
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Figure 5*~ Continued- 


MCA ]M L56G18 





(d) 011 = 20°. 

Figure 5.- Continued. 
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(d) 9]j » 20°. Continued. 

Figure 5*“ Continued. 
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(d) 0|j - 20°. Continued. 


Figure 5*- Continued. 
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("b) Effect of 7j. 



(c) Effect of 0j(f. 

Figure 6.- Example of the effects of jet Mach number, ratio of specific 
heats of the jet, and nozzle divergence angle upon the shape of the 
jet boundary. Pj/P„ = 5 . 
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Mj = 1.5. 


resBxire ratio upon the initial inclination 
oundary for 0^. = 0°. 

















Flgtcre 8.- Exan^les of Jet 
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(a) Ibctrapolation to 9 n = 0 °. 
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("b) Coinparison of extrapolated and calculated iDoundary. 

Figure 9 .- Example of extrapolation to 0if = 0°. Mj = 2.0; 7 j = 1.400 

and Pj/Peo = 10- 
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